Probing Non-Gaussianity In The Cosmic Microwave Background Anisotropics: One 

Point Distribution Function 
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We analyze WMAP 3 year data using the one-point distribution functions to probe the non- 
Gaussianity in the Cosmic Microwave Background (CMB) Anisotropy data. Computer simulations 
are performed to determine the uncertainties of the results. We report the non-Gaussianity param- 
eter /nl is constrained to 26 < /nl < 82 for Q-band, 12 < /nl < 67 for V-band, 7 < /nl < 64 for 
W-band and 23 < /nl < 75 for Q+V+W combined data at 95% confidence level (CL). 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Non-Gaussianity is one of the most important tests of 
models of the inflation. Among the various theoretical 
models on the inflation, slow-roll inflation is currently 
most lively being studied. There are various predictions 
on the magnitude of non-Gaussianity based on the sim- 
ple model of slow-roll inflation and its extensions, ranging 
from undetectably tiny values to large enough values to 
be detectable with currently available data [3, 0, 0, 0, S] • 
On the other hand, observational works have claimed 
both detection and non-detection of non-Gaussianity (for 
reviews on recent works, see H, 0, I, I E3)- Among 
the popular techniques for detecting non-Gaussianity 
are one-point distribution function fitting, bispectrum, 
trispectrum and Minkowski functionals. Here, we inves- 
tigate the one-point distribution functions to probe pri- 
mordial non-Gaussianity in the CMB anisotropy data. 
An observed CMB anisotropy at a direction (ST b s ) can 
be regarded as the superposition of three parts: physi- 
cal fluctuation of cosmic origin (5T P ), instrumental noise 
(<5T„), and foreground emissions (T/ ff ). Since the fore- 
ground templates are separately prepared, we start with 
foreground-removed data of which the CMB anisotropy 
can be decomposed into two uncorrelated components, 



ST = 6T obs - Tfg = 8Tp + ST, 



(1) 



The primary source for the cosmic fluctuation of CMB 
at the large scale is attributed to the Sachs- Wolfe ef- 
fect which is again triggered by the primordial curva- 
ture perturbation. The curvature perturbation $ by pri- 
mordial seed during the inflation is transferred to CMB 
anisotropy with the relation 



ST p (x) 
To 



r/t$(x) 



(2) 



where To = 2.725 K, the thermodynamic temperature of 
the CMB today, and r\ t is the radiation transfer function. 
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For the super-horizon scale, we take rjt = —1/3 from the 
Sachs- Wolfe effects. At the first-order of perturbation, we 
may replace $ — $ s , where $ s is an auxiliary random 
Gaussian field with its mean (& g ) = and its variance 
denoted by ($g)- When the second-order perturbation is 
considered, it is conventional to prescribe the nonlinear 
coupling of the curvature perturbation as [llT | 



$(x) 



$ 9 (x) + / NL ($»(x) -($»)) 



(3) 



where /nl is the non-Gaussianity parameter. The second 
term in ([3]) is responsible for the non-Gaussianity of the 
primordial fluctuation. Then, the probability distribu- 
tion function of the non-Gaussian field $ can be derived 



as 



/•(*) = / f G ^ g )S D [*-*„- /nl (*2 - d ®s 



nl(*+ 



exp 



2(1> 2 9 ) 



exp 



$2 



2<* 2 > 



where $± are defined by 



nl 



-l±,/l + 4/ NL $ + 4/ NL ($2) 



and $ has to be limited by the reality of <I>± as 



f * 1 
Jnl$ > -- 



($ 2 ) can be expressed in terms of rjt, Tq and acMB, 



XL 




(4) 



(5) 



(6) 



(7) 



For a pixelized CMB anisotropy data set, the probabil- 
ity distribution function for Gaussian instrumental noise 
becomes 



fN(o~T n 



1 



N„ 



A 



E 



i 



pix j~[ ^/2nal/n. 



exp 



6T? 



2er 2 M 



(8) 
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where rij is the effective number of measurements at the 
ith pixel and <To represents the dispersion of the instru- 
mental noise per observation (00 =2. 1898, 3.1249, 6.5112 
mK for Q, V, W-band, respectively Now, it is 

straightforward to express the probability density func- 
tion for ST in an integral form, 

f(ST) = J f S T p (ST p )f N (ST n ) 

xd D (ST - ST P - 5T n ) dST p d5T n 

= J M<S>)f N (6T n ) 

xd D (ST - r) t T <I> - 5T n ) d<5>dST„. (9) 

The probability density function derived in @ explic- 
itly contains the non-Gaussianity parameter /nl, and it 
can serve as the prediction of one-point distribution func- 
tion with a given /nl for a (ideally) foreground-removed 
CMB anisotropy data set to estimate the magnitude of 
deviation from Gaussian distribution in a quantitative 
manner. 



which is directly calculated from the WMAP data. Fig- 
ureCQand TableUshow the results of x 2 fitting of WMAP 
data sets with varying /nl as a free parameter. All data 
sets are best fitted at positive /nl (dubbed /nl ) which 
are consistent with one another as well as the results with 
previous work 0] within the statistical errors. 
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TABLE I: The results of x 2 -tests for Goodness-of-fit for 
WMAP 3 year data. The percentages on the first column rep- 
resent the tail probabilities at which the statistic \ 2 would be 
smaller than the observed. DOF on the bottom row stands 
for Degrees of freedom. 



II. APPLICATION TO WMAP DATA 

We use the three channels of WMAP 3 year CMB 
anisotropy data sets (Q (33GHz)-, V (61GHz)-, W 
(94GHz)-band) which contain dominant signal over con- 
taminations to investigate the non-Gaussianity of the 
CMB anisotropy data. To remove the foreground emis- 
sions, the Maximum Entropy Method (MEM) maps of 
the synchrotron, free- free and thermal dust are used[il|. 
The KpO-mask is applied to the sky maps to remove 
the intense Galactic emissions and scattered bright point 
sources, which leaves 76.5% of the sky. We also prepare 
a combined map (Q+V+W) by taking a weighted sum 
for a pixel temperature, 



5T(x) 



Ei n i ( x ) /°oi 



V, W (10) 



where rij (x) is the effective number of measurements at 
the pixelized position x and ao% is the dispersion of the 
instrumental noise of the ith channel. We can trace the 
effective variance of the instrumental noise as a result of 
weighted sum defined in (fT7j)) as 



(x) 



(x)/4 



Q, v, w. (li) 



The sky map data are degraded from N s idc = 512 (6.87' 
pixel) to N s ide = 128 (27.48' pixel) where the number 
of pixels in a full sky map is given by 12 x A 2 dc . The 
purpose of demotion of the resolution is to suppress the 
small scale fluctuation which is dominated by the instru- 
mental noise. We perform the x 2 -test for the goodness 
of fit for the probability density function given in §§§ as 
a prediction to the observed probability density function 



III. SIMULATION AND STATISTICAL 
UNCERTAINTIES 

As is shown in Figure [TJ WMAP data fit well with fi- 
nite range of the non-Gaussianity parameter /nl- We 
pick /^° L pt ' as the representative magnitude of non- 
Gaussianity for a data set, and carry out the Monte- 
Carlo simulations to test the pertinence of /nl as a 
proper measure of non-Gaussianity for a data set in a 
quantitative manner. A simulated data set is prepared 
as follows: First, a Gaussian field $ g with its variance 
equal to ([7]) is generated and we use it to generate 
field of which the deviation from Gaussianity is denoted 
by /nl- Second, we prepare a noise map in which each 
pixel contains a random value picked from a normal dis- 
tribution with variance Cg/n^ as defined in ||SJ) and add 
this to <I>-field. A simulated map generated in this way 
has the same noise structure and the dispersion of phys- 
ical fluctuation (ctcmb) as a real data set. Since a data 
set contains nontrivial instrumental noise and the num- 
ber of pixels is finite, the returned magnitude of non- 
Gaussianity (/nl^ ) would have some uncertainty. For a 
given value of /nl, we repeat the simulation and find that 
the algorithm returned an unbiased, normal distribution 
of /nl'*'' which is centered at the input value of /nl- 
Thus, from the results of the simulations, we are able to 
set the bounds on /nl for the real data. The simulation 
results and error bands are plotted in Figure and the 
deduced uncertainties for /nl are summarized in Tabic 
HT1 It is very remarkable that this analysis strongly disfa- 
vors the null hypothesis (/nl = 0) and all the data sets 
show consistent results within the statistical errors. The 
results are also consistent with the works by the WMAP 
team (-54 < / NL < 114 at 95% CL from bispectrum Q) 
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FIG. 1: x -test for goodness of fit is used to find the optimal value of the non-Gaussianity parameter. The anisotropy maps 
are demoted to iV s ide = 2 7 . Top left: Q-band, top right: V-band, bottom left: W-band, bottom right: Q+V+W combined. 
The color (shaded) bands indicate the bounds on /nl in which the probability that the statistic \ 2 is smaller than that of the 
observed with respect to the prediction curve with /nl. 



but with much tighter limits and more importantly, it 
excludes /nl = at 95% CL. 
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'Number of simulations for each input /nl- 

TABLE II: Summary of results from Simulations with WMAP 
3 year data profiles and the bounds at three confidence levels 
from the simulations. 



IV. CONCLUSION 

We developed an algorithm that uses the one-point dis- 
tribution function to investigate the non-Gaussianity of 
CMB anisotropy data, and applied it to WMAP 3 year 
data. We found that the null result (/nl=0) is mani- 
festly excluded at 95% CL. The estimated magnitude of 
non-Gaussianity parameter is 23 < /nl < 75 at 95% 
CL and 9 < /nl < 88 at 99% CL for the (Q+V+W)- 
combined map. Since the quadratic term in ^ takes 



a generic form of Taylor series for a perturbative ex- 
pansion, it is a good possibility that the observed non- 
Gaussianity in this work is a combined effects of various 
physical processes, while the primordial seeds are very 
likely to be the leading one. There are two premises we 
have taken in developing the algorithm, which, provided 
they are not precise enough, could cause non-Gaussianity 
of not cosmic but systematic origin: (1) the probability 
distribution function of the instrumental noise for each 
pixel is centered at zero, and (2) the foreground emis- 
sions are removed efficiently enough in the foreground- 
removed maps. The first condition can be broken when 
the thermal and radiation environments of the WMAP 
satellite in its orbit are taken into account, while the 
WMAP team assessed they are insufficient to influence 
the science data [La |. So, we tested the effects of the al- 
ternative noise distributions with a random mean in each 
of the Gaussian distribution in ([5]) and the algorithm was 
not misled to show non-Gaussianity within the statisti- 
cal error. It is difficult to directly estimate how much 
residual foreground emissions after foreground subtrac- 
tion would affect the one-point distribution function. We 
solely rely on the quality of foreground templates and it 
is remarkably successful, showing that the observed to- 
tal Galactic emission matches the model to less than 1% 
[HI, 03]. We also analyzed simulated maps which are 
(Gaussian map + foreground templates), and all the tern- 
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FIG. 2: Error bands found through the simulations with the same profiles with WMAP 3 year data (blue-green slanted bands). 
Here the "same profile" means that the simulated data set has exactly same instrumental noise structure and cfcmb as a real 
WMAP data set has. The vertical bands indicate the bounds for found from the analysis of the WMAP data. Top left: 

Q-band, top right: V-band, bottom left: W-band, bottom right: Q+V+W. On the vertical axis, /^J™'^ is the measured 
from the simulated data. 



plates for Q, V and W-channel showed negative values of 
the non-Gaussianity parameter with |/nl| ~ O (10 1 ) at 
the resolution iV s id = 512. 
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